NASA TECHNICAL 

MEMORANDUM 


so 

CO 

Cs| 

i 

X 


un 



LOAN COPY: RETUTF£ jg 
AFWL (DOUL) ^ = 

KIRTLAND AFB.tM^S 


COMPUTER METHOD FOR IDENTIFICATION 
OF BOILER TRANSFER FUNCTIONS 

by Jeffrey H. Miles 

Letvis Research Center 
Cleveland, Ohio 44135 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • DECEMBER 1971 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 


111 III HIM 


1. Report No. 

N ASA TM X -2436 

4. Title and Subtitle 


2. Government Accession No. 


nisi' 


b? 


5. Report Date 


COMPUTER METHOD FOR IDENTIFICATION OF BOILER 
TRANSFER FUNCTIONS 

7. Author(s) 

Jeffrey H. Miles 


9. Performing Organization Name and Address 


December 1971 

6. Performing Organization Code 


8. Performing Organization Report No. 

E-6485 

10. Work Unit No. 

112-27 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

12. Sponsoring Agency Name and Address 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


National Aeronautics and Space Administration 
Washington, D. C. 20546 

15. Supplementary Notes 


14. Sponsoring Agency Code 


16. Abstract 

An iterative computer procedure is described for identifying transfer functions using frequency 
response data. This procedure is applied to the study of boiler dynamics. The method uses 
the frequency response data to obtain a satisfactory transfer function for both high and low 
vapor exit quality data. 


17. Key Words (Suggested by Author(s)) Identification; Com- 
puter method; Frequency response; Transfer 
functions; Dynamic models; Boilers; Flow sta- 
bility; Two phase flow; Fluid dynamics; Inlet 
flow; Mathematical models; Instability 

18. Distribution Statement 

Unclassified - unlimited 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of Pages 

31 

22. Price* 

$3.00 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 



COMPUTER METHOD FOR IDENTIFICATION OF BOILER TRANSFER FUNCTIONS 

by Jeffrey H. Miles 
Lewis Research Center 

SUMMARY 

An iterative computer method is described for identifying boiler transfer functions 
using frequency response data. An objective penalized performance measure and a 
nonlinear minimization technique are used to cause the locus of points generated by a 
transfer function to resemble the locus of points obtained from frequency response 
measurements. Different transfer functions can be tried until a satisfactory empirical 
transfer function of the system is found. 

To illustrate the method, some examples and some results from a study of a set of 
data consisting of measurements of the inlet impedance of a single tube forced flow 
boiler with inserts are given. 


INTRODUCTION 

Studies of the dynamic stability of single tube boilers have been conducted at the 
NASA Lewis Research Center. The aim of these studies was to establish dynamic sys- 
tem stability criteria. The experimental part of the program consisted of measuring 
the frequency response of the system (refs. 1 to 3). As part of the analytical program 
a method was developed to identify transfer functions from these measurements. This 
report presents the details of this method. Typical results obtained by applying the 
method to the boiler frequency response data are also presented. 

The frequency response data consist of the amplitude ratio and the phase angle dif- 
ference of two parameters at various frequencies. The frequency response is con- 
sidered to be identified if a transfer function that fits the frequency response data is 
obtained. Various approaches can be used to analyze these data to obtain a satisfactory 
transfer function. 

If nothing were known about the dynamic system, a heavily empirical approach 
could be adopted. For example, it could be assumed that the system transfer function 
is composed of pure time delays e"^i, gains kj, and associated time constants t^. It 



could further be assumed that the transfer function Z had the general form 



i=l 


where s = jo>. This curve fitting approach does not have generality. The constants are 
not related by theory to the steady- state condition and must be calculated for each test 
measurement. 

Another approach to the problem would be to attempt a rigorous solution of the 
partial differential equations governing the boiler. This approach assumes all impor- 
tant processes are known. 

The approach taken herein was to develop a method that could yield an analytical 
model. The method is used to obtain a transfer function with the following characteris- 
tics: The transfer function should describe the measured data; the transfer function 
should preferably be derivable from a physical theory; and the dimensionless gains, 
time delays, and time constants in the transfer function should be derivable from 
steady- state measurements or physical properties. This report is concerned with the 
development of the method and not with relating results to fundamental physical prin- 
ciples. 

In the next section the general method is explained. In the then following section, 
the results of applying the procedure to boiler inlet impedance measurements are given. 


METHOD 

The success of the method is determined partially by the extent the locus of points 
generated by the transfer function (LPG) resembles the locus of points obtained from 
frequency response measurements (LPM). The locus of points in either case is ob- 
tained by connecting the frequency response points in the order of increasing frequency. 
While it is necessary that this resemblance criterion be satisfied, resemblance is not 
a directly measurable quantity. A secondary criterion called an objective performance 
measure is used to provide a quantity that is directly measurable. 

A description of the manner in which this performance measure is used to investi- 
gate resemblance between LPM and LPG follows. Each physical description (based on 
experimental results or fundamental physical principles) thought to describe the funda- 
mental dynamic phenomena is used to produce a system transfer function corresponding 
to the one measured. 
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Each physical description yields estimates of the gains, time delays, or time con- 
stants in the transfer function. It also yields information on the possible ranges for 
these constants. For each transfer function the constants are selected to minimize the 
positive, quadratic performance measure. The form of the performance measure used 
also ensures that all the constants will be within their proper limits. In this manner 
each analytical model can be compared to measured frequency response data. 

The objective performance measure indicates which physical description is most 
applicable. The best LPG can be checked for resemblance by plotting and comparing 
with the LPM. 

In this manner a transfer function is constructed by an iterative procedure. At 
each step the results can suggest modifications in the analyses that are necessary for 
improvement. This procedure can be used to pinpoint the inadequacies in an analysis. 

The basic performance measure S is 

N r ,, 

s =EF s i’ (x »-wM 

i=l 


where s A = j coj and the { x } are the transfer function constants. The performance 
measure is the sum of the square of the magnitude of the separation of the measured 
impedance Z Qbg and the calculated impedance Z at all frequencies where measure- 
ments are made. Previous reports (refs. 4 to 6) have suggested that a determination 
of the approximate system transfer functions and parameters can be made by finding 
the transfer function constants which minimize the performance measure. Standard 
computer techniques, to be discussed later, are used to find the constants { x } which 
minimize the performance measure. 

The actual performance measure used herein was a modified form of the previously 
given basic definition. First, the basic performance measure S was changed to a 
weighted performance measure S w> This was done for various reasons. The weighted 
performance measure provides a simple means to obtain a performance measure that 
applies over a limited frequency range. As an example, if measurements were made 
over a range from 0. 001 to 4. 0 hertz and if only the range from 0. 001 to 1 hertz was to 
be covered by the transfer function, the weighting factors for frequencies larger than 
1 hertz could be set to zero. The weighted performance measure furthermore provides 
a simple way to minimize the effects of noise or of errors on the identification proce- 
dure. It also provides proper scaling. The weighted performance measure has the 
form 
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S w (k,m) 


N 

i=l 


Z k (Sj, { x(m ) } ) - Z m (s.) 


The k index represents the theory. The m index represents the rrP 1 set of 
data. The weighting factor is designated as a(k, m)^. 

Since the selection of values of the transfer function constants is automated, one 
further modification was made to the performance measure. This modification was 
made to insure that the constants selected in minimizing the performance measure are 
in the proper range. The estimates of the proper range for each constant are used to 
provide constraints on the selection of parameters, and an acceptable variation of 
parameters is obtained. Hence, a transfer function need not be discarded if the initial 
selection of parameters was incorrect. An inadequate understanding of the process or 
a lack of knowledge of important steady- state parameters need not prevent the recogni- 
tion of a suitable transfer function if the range of the parameters can be estimated. 

The selection of parameters was constrained by the penalty function 


P ^ (x) = jd(k,m) z - x(k,m)J 2 S^ 

where 

fO for x(k,m)^ <d(k,m)^ 

= | 

Ij. for x(k, m)^ > d(k, m)^ 

This penalty function has the following characteristics. If the estimated parameter 
th th 

x(k, m)^ for the k Ln theory and in steady- state condition is less than the estimated 

constraint, the penalty p^ is zero. If the estimated parameter x(k, m) is greater than 
the estimated constraint, the penalty is equal to the square of the difference between the 
penalty and constraint. The penalty due to these individual penalty functions was defined 
as the weighted sum of these penalty functions. 

With this penalty function, the objective performance measure S was defined to 

r 

be 


N 

S p (k, m) = S w (k, m) + w^x) 

i=l 
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The weighting factor Wj is selected to scale the penalty such that each penalty will 
have an appropriate effect on S p . 

The problem of minimizing a nonlinear function like the performance measure can 
be solved by a search technique. The basic methods available are discussed in chap- 
ter 6 of reference 7. To provide maximum flexibility in the choice of Z, a search 
technique which does not require evaluation of derivatives was chosen. The search 
technique used was that of Powell (ref. 8). The computer program used was adapted 
from reference 9. 

When the penalized performance measure and the computer methods mentioned pre- 
viously are combined, the computer method shown schematically in figure 1(a) results. 
This method is used to find out if a particular choice of transfer function fits the meas- 
ured frequency response data. 

The procedure used to obtain a satisfactory transfer function is shown in figure 1(b). 
Only the step in which S (k, m) is calculated is programmed on the computer. For the 
first transfer function it is best to select one based on a simple first-order theory that 
yields an approximation to the LCPM and that uses the steady-state data available to 



(a) Procedure used to check transfer function Z^s, x). 
Figure 1. - Flowcharts illustrating identification procedures. 
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(b) Procedure used to find a satisfactory transfer function. 
Figure 1. - Concluded. 


predict the transfer function constants and constraints. At this point it is best to select 
the data that are easiest to work with and to limit the frequency range over which agree- 
ment is expected. This initial transfer is denoted as Z^. The set of parameters ob- 
tained is denoted as { x(l, m) } and the set of constraints is { d(l, m) } . The minimum 
value of Sp(l, m) can then be found. It should be noted that in this expression Sp may 
be minimized to yield a transfer function Z j which is optimum over a limited fre- 
quency range. 

For the following iterative steps it is necessary to examine previous results to 
locate terms that would extend the frequency range, terms that would enable the theory 
to be extended to more steady- state operating conditions, and terms that were not used 
because of insufficient steady-state knowledge or physical knowledge. After the func- 
tional form of these neglected terms is found, the value of the parameters in these 
terms and the constraints on these parameters are then estimated. These terms when 
combined together produce various theoretical transfer functions. The k Ln theoretical 
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description yields transfer function Z^. Also from the k theoretical description and 
for each steady- state condition, a set of constants { x(k, m) } and a set of constraints 
{ d(k,m) } are available. The minimum value of S p (k,m) is found. The resulting set 
of optimized transfer functions Z* is compared. 

The results of each change in the theoretical transfer function, change in frequency 
range, or change in constraints on the previous theoretical transfer function are readily 
available using the computer method. 

The procedure described previously was used to obtain a transfer function for the 
inlet impedance of a single tube, forced flow boiler with inserts from inlet impedance 
measurements. The following section discusses some of the results of this work. The 
aspects discussed in the next section are those not concerned with relating the transfer 
function obtained to a set of fundamental physical principles. 


APPLICATION 

This section initially indicates how the contributions of various transfer functions 
(each attributed to a different physical process) are separated in the formulation of the 
system transfer function. This is done so that each can be considered independently. 
The transfer function for one process can then be changed until the system transfer 
function is improved while maintaining the transfer functions for other processes con- 
stant. Next, some typical results obtained by using the performance measure, com- 
puter techniques, and the iterative procedure described in the previous section are 
given. Last, the transfer functions which were found by this method are discussed. 

The analyzed data (ref. 2) were obtained experimentally by driving sinusoidally the 
open area of a valve in the feed system about a mean area and measuring the pressure 
and flow at the boiler inlet. The pressure and flow signals were analyzed by a 
frequency- response analyzer. This analyzer computed the magnitude and phase of the 
sinusoidal content of the pressure and flow signals at the driving frequency. The mag- 
nitude of the boiler inlet impedance at each frequency is equal to the amplitude ratio of 
the inlet pressure to inlet flow signals. The phase of the boiler inlet impedance at each 
frequency is equal to the difference between the measured phase angle of inlet pressure 
and inlet flow. 

The theoretical model used to describe the measured boiler input impedance data is 
a function of a combination of physical processes. The different physical processes are 
assumed to be tied together by a set of equations. While the detailed models of each 
physical process are not known adequately, the set of equations which forms the frame- 
work that ties the models together is assumed to be known. 

The equations and subsystem transfer functions are derived from the governing 
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equations for the boiler by (1) simplifying and linearizing for small perturbations about 
a known steady- state operating condition, (2) Laplace transforming the equations, and 
(3) putting the equations in dimensionless form. The dimensionless form is obtained by 
dividing by scale factors P and W. The scale factors are not defined individually 
since in the final equation they appear as W/P. The quantity P/W is defined to be 
3.28x10° newtons- seconds per square meter - kilogram. The result is a set of equa- 
tions connecting perturbations in parameters. The coefficients of the parameters are 
the subsystem transfer functions. Each subsystem transfer function represents a 
physical process. 

The perturbation parameters chosen to describe the inlet impedance are the follow- 
ing: the variation of inlet pressure AP. , the variation of inlet mass flow rate AW. , 

o * 111 in ? 

the variation of vapor pressure beyond the subcooled region APx , the variation of 

sc 

subcooled length AL and the variation in vapor mass flow rate beyond the subcooled 

b C 

region AWt 

^sc 

These five perturbations are connected by the following four equations: 


AP. 


AW- AP L 

™ = E (s)-^ + -_- s .c 


W 


AL AW. AP L 

^gg = G(s)^g + H (s) - -, SC 

L sc W P 


( 1 ) 

(2) 


AL 


AW L AW. 

W W * L 


AP, 


sc 


F(s) • 


sc 


SC 


(3) 


ap l aw l 

_SC=R __SC (4) 

P W 

Equation (1) relates the perturbation in pressure drop across the subcooled region to the 
inlet mass flow rate through an inertial term E(s). Equation (2) states the following: 

(1) Changes in subcooled length (for constant vapor pressure) are due to perturba- 
tions in inlet mass flow rate effecting heat transfer (the transfer function for this phys- 
ical process is represented by the term G(s)). 

(2) For constant inlet mass flow, perturbations in vapor pressure effect the sub- 
cooled length (the transfer function for this physical process is represented by the term 
H(s)). 

(3) These two effects can be added for small perturbations. 
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Equation (3) states that changes in mass flow rate can be related to mass storage 
through changes in subcooled length (represented by the term figs) and in pressure 
(represented by the term F(s)). Equation (4) relates variations in exit flow and exit 
pressure by the constant R. 

The inlet impedance derived using these equations is 


Z(s) = 


AP m/ P R ( 1 - a 2 s °( s >) 

A w in /W [l + R { 0 2 sH(s) + F(s) }" 


+ E(s) 


(5) 


The basic block diagram for the system is shown in figure 2. 

Each subsystem transfer function and thus each separate physical process can be 
investigated independently by proceeding in the previously stated manner. Thus, differ- 
ent theoretical models for each physical process can be checked. Also, different math- 
ematical forms of the subsystem transfer functions can be tested. These forms can be 
those determined from physical theory or from transfer function plots (ref. 10). The 
method of this report was then used to determine the subsystem transfer functions G(s), 
E(s), H(s), and F(s). In the following section the results of a study of data taken at a 
condition of high vapor exit quality are discussed. 


High Quality Vapor Exit Condition 

The case discussed is run 4 from reference 2. The data were obtained at a steady- 
state condition of high vapor exit quality (99 percent). The boiler inlet impedance data 
are plotted in figures 3(a) and (b). Figure 3(a) contains the data measured over the 
frequency range 0. 04 to 1.0 hertz, and figure 3(b) contains plots of the data over the 
frequency range 0. 25 to 4. 0 hertz. The basic characteristic form is that of a spiral. 
Over the frequency range below 1 hertz (fig. 3(a)) the amplitude of the transfer function 
decreases as the phase angle rotates clockwise. Above 1 hertz (fig. 3(b)) three impor- 
tant changes in the form occur. Between 1. 6 and 2. 5 hertz a small loop occurs. Be- 
yond 2. 5 hertz after the small loop the following loop is larger. Also above 4 hertz the 
data are not symmetric about (0, 0). Instead, the data are symmetric about a line 
through the real axis. Also, the spiral seems to be pulled up along a line parallel to the 
imaginary axis. Several theories are now compared with the data. 

Theory 1 . - The initial transfer function used to provide the starting point in the 
iterative model building process was the theoretical model described in reference 11. 
This model is designated as theory 1. This theory, which is based on many simplifying 
assumptions, provides an adequate estimate of the boiler impedance over a limited fre- 
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queney range for those steady- state conditions where the theoretical simplifications 
agree with the actual steady- state conditions. The theory was not meant to be appli- 
cable to data taken at a condition of very low exit quality or very high exit quality. 

The model described in reference 11 used the following functional forms for the 
subsystem transfer functions G(s), E(s), H(s), and F(s): 

E(s)=x 4 s (6a) 

-& 2 S 

G(s) (6b) 

e 2 s 

H(s) = x 2 (6c) 

F(s) - 0 (6d) 


The time constant 0 2 


is the subcooled dead time and is given by 



L scP A 

w 


(7) 


The time constant 0 2 is designated as constant x^. The time constant x 4 corresponds 
to the inertia of the subcooled fluid. It is given in the dimensionless form used here as 


x 4 ~ 


L sc W 
A P 


( 8 ) 


The gain Xg corresponds to the change in subcooled length with a change in vapor pres- 
sure. In nondimensional form it is 


0L 


x 2 = 


sc 


9Pt 


sc 


sc 


(9) 


The boiler resistance previously denoted R is henceforth defined as a constant: 


x 3 =R 


(10) 
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To obtain a form of the inlet impedance equation that is more clearly related to more 
common forms, the following transformation of equation (5) was made. Using equa- 
tion (10) and defining 


QUA = 1 - Xl sG(s) 


( 11 ) 


give 


XoQUA 

Z= *— — + E(s) (12) 

1 + XgXjSH(s) + XgF(s) 

The first transfer function used in the procedure diagrammed in figure 1(b) was 
derived from equations (11) and (12) and equations 6(a) to (d). The result is 


Z - 


x 3 e 


•XjS 


1 + XjjXgXjS 


(13) 


Note that from equations (11) and (6b) 


QUA X - e 1 (14) 

Coefficients for this model were selected using the method described in the first 
section. The resulting function is plotted for frequencies less than 0. 3 hertz in fig- 
ure 3(c) and for frequencies greater than 0.3 hertz in figure 3(d). Examination of fig- 
ure 3(c) shows that this model does describe the data in the range below 1 hertz. The 
data and the theory 1 locus both resemble a spiral. The amplitude of the transfer func- 
tion decreases as the phase angle rotates clockwise. Examination of figure 3(d) shows 
that theory 1 fails in the frequency range above 1 hertz. The sizes of the loops calcu- 
lated decrease uniformly while the measured transfer function shows a small loop fol- 
lowed by a large one. Also in the same figure the calculated spiral is seen to by sym- 
metric about the imaginary axis whereas the measured transfer function is seen to be 
symmetric about a line parallel to the imaginary axis on the positive real side. An 
aspect of the theory 1 curve and the data that are in agreement is the appearance of 
loops. 

The relation is now shown between equation (13) and the theory 1 curve in figures 
3(c) and (d). The spiral shown in figure 3(c) is generated by the term 
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( 15 ) 


x 3 e 


-XjS 


1 + XgXgXjS 


The decrease in amplitude occurs because the magnitude of 


1 + XgXgXjS 


(16) 


decreases with frequency. 

The change in phase angle with frequency is seen as a rotation in figures 3(c) and 
(d). The rotation is due to a slow change in phase angle from equation (16) combined 
with a rapid change in phase angle due to the term 

e " X l S (17) 


The addition of the imaginary term 


x ,s 


(18) 


to term (15) causes the spiral to be pulled up along the imaginary axis. This effect 
produces the loops seen in the theory 1 curve of figure 3(d). Theory 1 can be modified. 
The next sections discuss one type of modification. 

Theory 2 . - As an example of one step in the iteration procedure, the results ob- 
tained by modifying the previous model to include the effect of wall storage on the heat 
transfer process are discussed. The wall storage modification is typical of a modifica- 
tion that can usually be made from a more detailed analysis of the physical processes 
in a system. 

The modification results in a change in the definition of QUA. This modification 
will be called theory 2. The modified term was 

QUA 2 = ;ue" X + 1 - ju (19) 


where 


l-t — 1 + 


x 1 s(l - x 13 /x 1;l ) 

1 + x 12 /x 11 + XjS/Xjj 


(20) 
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( 21 ) 


X-.S + X 


X =- 


12 " x 12 


1 + XjS/x^ 


Thus, using equations (12), (6a), (6d), and (19) yields 


'1A 


(jne - ^ + 1 - p)x, 

1 + XjXgXgS 


+ x 4 s 


( 22 ) 


A range of parameter values for x^, x-^j and x^g was selected. Also, values of 
the coefficients for the specific operating condition of the example were chosen. The 
resulting transfer function is shown in figures 3(e) and (f). 

Examination of figure 3(e) shows that the theory 2 curve describes the data at low 
frequencies as well as the theory 1 curve. The procedure used has minimized the per- 
formance index. Still, a glance at figure 3(f) indicates that the theory is not successful 
since the resulting curve does not resemble the data at high frequencies. The theory 2 
loops decrease in size rapidly as the frequency increases. Also, no small loop followed 
by a large one occurs. This modification is an example of how the extension along well 
known lines of an idea may result in an unsuccessful transfer function. The final ver- 
sion is discussed next. 

Theory 3 . - The model found to represent the high vapor exit quality data best con- 
sisted of the following set of subsystem transfer functions: 


E 


QUA = e 


-XjS 


= x 4 s + 


l/x iQ + S/Xg 


H = XgXjSe 


-x 5 s 


F =- 


X 8 S 


S + XgS + Xrj 


(23) 

(24) 


(25) 


(26) 


The resulting transfer function, which is plotted in figures 3 (g) and (h), is desig- 
nated as theory 3. This model has all the important characteristics of the data at all 
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frequencies. The transfer function has the basic spiral shape below 1 hertz (fig. 3(g)). 
It also has the small loop between 1. 6 and 2. 5 hertz. And the larger loop between 2. 5 
and 4 hertz also exists (fig. 3(h)). Furthermore, the high frequency data are symme- 
tric about a line through the real axis. The inertia effect that pulls the data up parallel 
to the imaginary axis is also apparent. 

The transfer function is 


-XjS 


J 1B 


-XgS 

1 + x-^xgse + 


XgXgS 


+ x 4 s + 


l/x i0 + S/Xg 


2 2 
S + XgS + Xrj 


(27) 


The appearance of the small loop before the large one (fig. 3(d)) is due to the addi- 
tional term F(s). The magnitude of the term 


-x 5 s 

1 + XjXgXgSe + 


X 3 X 8 S 


2 2 
S + XgS + Xrj 


(28) 


becomes very small at the frequency where 

s 2 + X 2=0 

This decrease in magnitude causes the loop that occurs between 1. 6 and 2. 5 hertz to be 
small. The high frequency data are symmetric about a line through the real axis paral- 
lel to the imaginary axis due to the term 


l/x 1Q + s/x 9 


(29) 


in E(s). At high frequencies including the term causes E(s) to become 

E(s) = x 4 s + Xg (30) 

In conclusion, it should be noted that before equation (27) was determined no equa- 
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tion existed which gave good agreement with the form of the high vapor exit quality data 
for the frequency range above 1 hertz. 

The three theoretical models used in the study of the high vapor exit quality example 
are summarized in table I. The set of coefficients that was obtained using the method 
with each model and the data are given in table II. Some additional results of the study 
are presented next. 


TABLE I. - COMPARISON OF FUNCTIONAL FORMS 

OF TRANSFER FUNCTION TERMS 

FOR DIFFERENT THEORIES 

|7i = 1 + XjS(l - x 13 /x n )/(l + x 12 /x 11 + x 1 s/x n ); 

X = XjS + Xjg ~ x^ 2 /(l + XjS/x^)T] 


Transfer 

Theory 

function 

term 

1 

2 

3 

QUA 

-x^s 

e 

Me"* + 1 - m 

-x-s 
e 1 

El 

X 4 S 

X 4 S 

X 4 S 4- s/(l/x 10 + s/Xg) 

HI 

X 2 X 1 S 

X 2 X 1 S 

-X 5 s 

x 2 XjSe 0 

F3 

0 

0 

x 8 s/(s 2 + 2sxg + x 2 ) 


TABLE n. - DIMENSIONLESS TRANSFER FUNCTION 
GAINS 3 AND TIME CONSTANT PARAMETERS 


Transfer function 
parameter 


Theory 

12 3 


x i> 

sec 

x 2 


x 3 


x 4 , 

sec 

x 5’ 

sec 

x 6 , 

sec" 

x 7 , 

sec" 

x 8> 

sec” 

x 9 


x id 

, sec 

X 11 


x 12 


x 13 



0.794 

1.96 

0.855 

l.OOxlO ' 2 


0.706 

1.72 

0.881 

5.00xl0~ 3 


190.0 

3.06 

1.46 


0.883 
0.44 
0.845 
0.00296 
0. 174 
2. 01 
12.8 
138.0 
0.0365 
0.098 


a Gains relative to 3.28x10^ N-sec/m 2 -kg. 
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Low Quality and Superheat Vapor Exit Condition 


No equation for the LPM of the low vapor exit quality data existed before it was 
determined that the following modification of equation (27) resulted in a good fit. Equa- 
tions (23) to (25) are still used in equation (12). But, instead of equation (26) the follow- 
ing form of F(s) is used: 


SXo 

F(s)= — (31) 

1 + SXj J 

The excellent agreement possible with this transfer function is demonstrated for a typi- 
cal low vapor exit quality (20 percent) case (ref. 2, case 1) in figure 4. 

When the method was applied to data taken at an exit condition of superheated vapor 
(ref. 2, case 9), it was found that equation (27) produced good results. A comparison 
of these results to the data is made in figure 5. Excellent agreement between the meas- 
ured and calculated phase and amplitude over the whole range of frequencies from 0. 01 
to 4 hertz exists. 

The procedure described in the report enabled a determination of the transfer func- 
tion form to be made for low exit vapor quality, high exit vapor quality, and superheated 
vapor exit conditions. Thus, the procedure devised has proved to be valuable in extend- 
ing and increasing the understanding of the dynamics of boilers. 


CONCLUDING REMARKS 

The search technique used permitted the parameters which were derived from the 
most certain theoretical aspects of the problem to stay relatively unchanged during the 
search for a minimum of S p . A good example was the term Xg which is calculated 
from the slope of the pressure drop against flow curve (ref. 11). This term determined 
the magnitude of the transfer function at low frequencies, and it varies by less than 
2 percent in a typical computer calculation. 

The parameters which were less certain took on a larger range of values during a 
search. The difficulty is that the Powell method generally provides only a local mini- 
mum and this means that the solution is not inherently unique. The theoretical results 
are still considered valid since they apply to a range of experiments with different 
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steady- state conditions and boiler geometries. The value of the constants obtained are 
also physically reasonable since the use of penalty functions provides estimates which 
are possible. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 14, 1971, 

112-27. 
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APPENDIX - SYMBOLS 


A 


d 


E(s) 

F(s) 

f 

G(s) 

H(s) 

I 



al/l 

sc' 


sc 


P 


P 

AP/P 


P 

R 

S 


S, 


v 


s 

T 

W 


cross-sectional area of 
2 

boiler tube, m 

constraint on value of trans- 
fer function constant 

subsystem transfer function 

subsystem transfer function 

frequency, Hz 

subsystem transfer function 

subsystem transfer function 

inertia, sec 

v=l 

gain 

subcooled length, m 

perturbation in subcooled 
length 

9 

mean pressure, N/m 

9 

pressure, N/m 

perturbation in pressure 

penalty function 

subsystem transfer function 

performance measure 

penalized performance 
measure 

weighted performance 
measure 

Laplace operator, sec - '*' 
time delay, sec 
mean mass flow, kg/sec 


W 

AW/W 

w 


{ x} 

Z 

Z = (P/W) 


a 


5 


mass flow, kg/sec 

perturbation in mass flow 

weighting factor in penalty 
function 

set of transfer function 
constants 

transfer function 

impedance value relative to 
which all transfer func- 
tions are compared, scale 
factor (3. 28x10® N-sec/ 
m 2 kg)" 1 

weighting factor in perform- 
ance measure 

weighting factor for penalty 
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(a) Boiler data for frequencies belov 1 hertz. 

Figure 3 ■ Comparison ol boiler data and eiperimenial models lor condition ol 99 percent vapor enl qualify fref. 2 run «). (Impedance 
relative to 3 28i Hr N - setlm‘-kg. i Note; To avoid confusion and increase graphical clarity, frequency calculatedirom the theories 
was nrt shown on the curves While judging agreement in magnitude is possible, judging agreement in phase angle is nol. To see 
that at an frequencies it is possible to achieve good agreement of Ihe phase angle and magnitude for the data and Ihe calculations 
based on theory 3. see figures e ana 5 Phase and magnitude for curves based on thrones I to 3 in (his figure agreed with the 
data at frequencies up to l her!/, beyond 1 hertz only theory 3 Shows good agreement 





240° 250° 260° 270° 

(b) Boiler data for frequencies above 0. 25 hertz. 
Figure 3. - Continued. 
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(e) Comparison of boiler data and theory 2 at frequencies below 0.3 hertz. 
Figure 3. - Continued. 
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(f) Comparison of data and theory 2 for frequences above 0. 3 hertz. 
Figure 3. - Continued. 
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Phase of boiler inlet impedance, deg 




Figure 4. - Comparison of experimental data and calculated boiler inlet impedance for exit condition of 20 percent 
vapor exit quality (ref. 2, run 1). (Impedance relative to 3. 2xI0 6 N-sec/m 2 ~kg. ) 
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